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SUMMARY 


The results of an approximate nonlinear finite-element analysis of a single 
lap joint are presented and compared with the results of a linear finite-elanent 
analysis, and the geometric nonlinear effects caused by the load- path eccentric- 
ity on the adhesive stress distributions are determined. The approximate non- 
linear finite-element solution is evaluated by ccmparison with a nonlinear solu- 
tion obtained from a finite-difference analysis, and results are compared with 
the classical approximate analysis by Goland and Reissner, An experimental pho- 
toelastic study of the moment distribution in the lap-joint adherend is used to 
evaluate the ade^acy of plate approximations for analyzing single lap joints. 
The results from finite-element, Goland-Reissner, and photoelastic analyses show 
that for a single lap joint the effect of the gecmietric nonlinear behavior of 
the joint has a sizable effect on the stresses in the adhesive. The Goland- 
Reissner analysis is sufficiently accurate in the prediction of stresses along 
the midsurface of the adhesive bond to be used for qucditative evaluation of the 
influence of geometric or material parametric variations. Detailed stress dis- 
tributions in both the adherend and adhesive obtained from the finite-element 
analysis are presented to provide a basis for ccmparison with other solution 
techniques. 


INTRODUCTION 


Adhesively bonded joints offer the aerospace designer an attractive mass- 
efficient alternative to the mechanical fastening of structural components, 

(See, for example, ref. 1.) In the past, however, the poor reliability of 
bonded joints has in general dissuaded designers from taking full advantage 
of the projected benefits of bonded structural designs. Thus, growth in the 
use of bonded systems in load-carrying components is not keeping pace with the 
continuing development arei improved reliability of new high-strength adhesives 
(refs. 1 and 2). Bonded joints should be designed to transfer load in shear 
with a minimiffn of peel across the bond line since adhesives are generally more 
efficient in supporting shear forces and perform poorly when sipporting peel- 
type forces. Design of bonded joints thus requires knowledge of the shear and 
peel stress distributions in the bond line as well as the shear strength of the 
adhesive. The single lap joint (fig. V) is the simplest joint which transfers 
in-plane load primarily by shear in the adhesive. In fact, the ASTM D 1002-72 
single-lap- joint test (ref. 3) is the roost popular experimental procedure for 
evaluating the shear strength of adhesives. 

Application of an in-plane tensile load to the single lap joint causes the 
adhe rends to deform due to eccentricity in ths load path. The deformations that 
result are of the order of the thickness of the adherend and an order of magni- 
tude larger than the thickness of the bond line. As loading progresses, the 
internal moment at the edge of the overlap region is reduced considerably due 
to the deformation of the adherends. Since the moment directly influences the 



shear and peel stress distributions in the adhesiver nonlinear solution proce- 
dures which reflect this reduction in mcment may be necessary for accurate 
determination of these stresses, 

A classic approximate analysis of single lap joints that takes into account 
the nonlinear action of the lap joint due to flexure of the adherends was pub- 
lished in 1 944 (ref. 4) and is hereafter referred to as the Goland-Reissner 
solution. This analysis established the basis for the study of the stress dis- 
tributions in single lap joints. The analysis was based on several assuraptlons, 
the most constraining being the assumption that the shear and peel stresses are 
both constant through the adhesive thickness. The assumption of constant shear 
stress leads to a nonvanishing adhesive shear stress at the free edge of the 
adhesive. The desire to remove these assumptions and to include other consider- 
ations, such as adhesive plasticity and adherend transverse shear deformation, 
prompted several researchers to refine the analysis (refs. 5 to 14). These 
refined analyses are based on fundamental plane strain elasticity approaches 
which require some approximations to either the governing equations or the 
solution techniques. 

By the midsixties, researchers were approaching the problem using finite- 
element analyses. Most of the investigations were linear and used mesh sizes 
that were too coarse to adequately predict the adhesive shear-stress behavior 
near the free edge (refs. 15 to 20). One linear investigation did use a fine 
mesh to determine the f ree-edge conditions in the lap region (ref. 21 ) . Several 
investigations were conducted using a nonlinear finite-element solution proce- 
dure to determine the stress distributions in the plastic range of the adhesive 
material (refs. 22 to 24) but did not consider the geonetric nonlinear behavior. 

The main purpose of this paper is to investigate the effect of geometric 
nonlinearity on the elastic stress distributions for single lap joints. First, 
the Goland-Reissner solution with slight modifications to the governing equa- 
tions is compared with the same solution with the nonlinear terms suppressed. 
Next, an approximate nonlinear finite-element solution is obtained and its 
accuracy confirmed by comparison with a solution from an established nonlinear- 
analysis computer program which uses a finite-difference solution technique. 
These solutions are also compared with corresponding linear solutions to indi- 
cate the influence of the nonlinear terms. The finite-element solution is com- 
pared with the solution obtained from the Goland-Reissner analytical approach 
to investigate the limitations of this latter approach. Finally, detailed 
distributions through the adherend and adhesive obtained from the nonlinear 
finite-element solution are presented to provide a detailed picture of the 
stress behavior in a single lap joint loaded in the elastic range. The 
stresses thus obtained are examined to evaluate the appropriateness of the 
single-lap- joint configuration as a test specimen for the determination of 
adhesive shear strength. 

Certain commercial materials are identified in this paper in order to 
specify adequately which materials were investigated in the research effort. 

In no case does such identification imply recommendation or endorsement of the 
product by NASA, nor does it imply that the materials are necessarily the only 
ones or the best ones available for the purpose. 
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SYMBOLS 


A extensional stiffness of adherend defined after equation (A3) 

a cross-sectional area of adherend 

c one-half overlap length 

D plate stiffness of adherend defined after equation (A2) 

E Young's modulus 

G shear modulus 

I centroidal moment of inertia of cross section of photoelastic 

adherend 

k,k' constants defined after equations (A8) and (All), respectively 

% length of adherend from point of load application to beginning of 

overlap 

M adherend moment resultant 

(Mq) 2 adherend mconent resultant at edge of overlap for stiff adhesive 
model based on linear analysis (see fig. 30) 

P applied load on photoelastic model 

ri ,r 2 constants defined after equation (All) 

T adherend load resultant, force per unit width 

t adherend thickness 

u adhesive displacement in direction normal to plane of adhesive 

u-] ,U 2 constants defined after equation (A1 5) 

V adherend transverse shear-force resultant 

w lateral displacement of adherend 

x,y lineal measure defined where used 

A constant defined after equation (All) 

6 constant defined after equation (A6) 

^ adherend strain in longitudinal direction 

5 ratio of adhesive to adherend thickness 
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T) adhesive thickness 

X coefficient defined after equation (All) 

V Poisson's ratio of adherend material. 

o adhesive peel stress 

CTjj normal stress in x- direct ion 

Oy norm^ stress in adherend in y-direction 

T adhesive shear stress 

<t> angular measure between line of axis of applied force and plane of 

bond 

Subscripts: 
a adhesive 

app applied 

av average 

L lower adherend 

max maximum 

U upper adherend 

o adherend stress resultcuit at edge of overlap 

Whenever terms in appendix A differ from corresponding terms in refer- 
ence 4, these terms are underlined with a dashed line, 

GOLAND-EEISSNER ANALYTICAL SOLUTION 
Lap-Joint Configuration 

The single-lap-joint configuration used throughout this study is based on 
the ASTM D 1002-72 test specimen used to obtain strength properties of adhesives 
in shear by tension loading. The dimer^ions of the joint are given in figure If 
the material properties assumed for this study are given in table I and are 
based on a titanium adherend and an aluminum-filled polyimide adhesive. The 
adherend thickness used in this study is slightly less than that suggested by 
ASTM D 1002-72. 
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Adhesive Stress Distribution 


The classical approximate analysis of reference 4 for a single lap joint 
with a flexible adhesive is repeated in appendix A with some slight modifica- 
tions to improve the consistency of the analysis. The shear distribution x 
and the peel stress (stress normal to the plane of the bond) a in the adhesive 
along the bond line are given by 


1 


^ av ^ L 


(1 + 3 k) 


sinh 6 


/S x\ 

cosh — + 3(1 - k + ?) 
\c / _ 


(A8) 


and 


^ app ' c ^ ^ ^ 2 


+ k'X sinh X s 


in xj sinh^ 


X x\ / Xx 1 

sini — i 

c / \c / 


^ \ 

+ 1 — + k'X cosh X cos X 1 ooshl — ) cos! ^ 

AC 


V 2 


\c 


c U 


(All) 


where all quantities are defined in appendix A. The constants k and k' are 
linearly related to the applied load. 

In the derivation of these equations, several assumptions are made. The 
adhesive shear and peel stresses are assumed to be constant through the adhesive 
thickness, and the adhesive direct stress parallel to the 'adhere nds is assumed 
negligible. The adherend and adhesive materials are assumed to be elastic and 
to behave linearly. The shear deformations in the adherends and the influence 
of the adhesive on the flexural stiffness of the joint are assumed negligible, 
and the adherends are assumed to deform as a plate in cylindrical bending. 

As a result of these approximations, an adhesive free-edge boundary condi- 
tion, T « 0, is unsatisfied. Concern with the nonvanishing shear stress at 
the edge of the adhesive is probably unjustified since real lap joints often 
have adhesive runout at the overlap edges. Although this adhesive runout is 
probably of lesser strength and stiffness than the adhesive in the overlap 
region, which is cured under pressure, it can support shear stress. An analysis 
of the shear stress distribution in an adhesive fillet at the edge of the lap 
is studied in reference 21. The authors assumed that the adhesive fillet 
material had the same properties as the adhesive in the bond line. 
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Mherend Mcroent Distribution 


The transverse deformations of the adhe rend are shown in figure 2 for two 
load levels. In a linear analysis, these deformations would be ignored in com- 
puting the moment distribution along the adhe rend. When accounted for in the 
calculation of the moment, these deformations can reduce the moment considerably 
from that calculated by assuming no relieving deformation, as shown in figure 3, 
The moment distribution is obtained by using a plate theory approximation, the 
accuracy of which is questionable in the vicinity of the overlap because of 
large transverse normal stresses and shear stresses which exist in this region. 
In order to investigate .the accuracy of this approximation, a photoelastic 
study of a series of lap joints was performed and is presented in appendix B. 

The photoelastic study showed that the approximation for the mcanent distribu- 
tion in the adherend is better for lap joints when the adhesive is flexible 
compared with the adherend than when the adhesive has the same stiffness as the 
adherend. The results of the study indicate that the adherend does have a 
stress distribution consistent with that of a cylindrically bent plate, even 
near the overlap region. 


Comparison of Linear and Nonlinear Stress Distributions 


The computed adhesive nondimens ional shear and peel stress distributions 
along the bond line are shown in figures 4 and 5, respectively, for applied 
loads in the working range for common adhesives and for the linear solution 
which gives stress concentrations independent of load level. The computed 

stress concentrations at the edge of the lap /- = 0^ are reduced for both the 


\° J 

shear and peel stresses as the applied load increases and are considerably 
below the stress concentrations predicted by linear solution. Differences in 
adhesive stresses are less severe at points in the interior of the lap region. 
The stress concentrations are plotted as a function of applied load in figure 6. 
The linear solution gives the maximum values of the stress concentrations and 
represents the limiting case for an infinitesimal applied load. The stress con- 
centrations canputed with the nonlinear terms retained drop appreciably as the 
applied load increases from zero and then drop more slowly as the load increases 
further. 


As the overlap length increases, the difference between the linear and non- 
linear analyses increases rapidly; and as shown in figures 7 and 8, the linear 
analysis predicts stress concentrations considerably higher than those predicted 
by the nonlinear analysis. Several authors have investigated the accuracy of 
the Goland-Reissner nonlinear analysis with finite- element studies (refs. 15, 

17, 21, and 22). In each instance the studies were made neglecting the non- 
linear effects. 


NUMERICAL ANALYSES 

As discussed in the previous section, the Goland-Reissner analysis makes 
several simplifying assumptions. Finite- element and finite-difference analysis 
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computer codes that can more realistically model the physical conditions are 
available and have been jsed to analyze the ASTM D 1002-72 single lap joint, 
this section^ a quasi— nonlinear finite-element solution of the lap joint is 
developed, coRipared with a nonlinear finite-difference solution to confirm its 
accuracy, and is then used to indicate the importance of the nonlinear influ- 
ence. All numerical calculations were conducted for a specimen with the dimen- 
sions shown in figure 1. 

Nonlinear Finite-Element Analysis 

The finite-element program SPAR described in reference 25 was used to 
analyze the single lap joint. SPAR is intended to be used primarily to perform 
structural analyses of linear problems. The major motivation for using the SPAR 
program was its capability to efficiently handle large problems with a minimum 
of execution time, central memory, and secondary data storage requirements. 

These characteristics are particularly advantageous for the reanalysis proce- 
dures re<^ired in the nonlinear-analysis approach. 

The finite-element mesh used to analyze the single lap joint is shown in 
figure 9(a). A total of 702 nodes and 624 triangular and quadrilateral membrane 
elements are used in the analysis. The element formulation was based on an 
assumed stress field with a minimum complementary energy method. (See ref. 26.) 
A similar analysis for a model with a smaller number of elements (496) gave 
results only slightly different from the present model, thus indicating that 
the modeling used is adequate. A very fine mesh was used in the critical areas 
where large stress gradients are expected, such as at the end of the adhesive 
and in the surrounding adherends. (See enlarged detail, fig. 9(a).) The 
adherend and adhesive thicknesses were each divided into five equal layers. 

In high stress gradient areas, the aspect ratios of the elements were always 
maintained between 1.0 and 3.0 to assure the accuracy of the analysis. 

The closed-form analysis showed the importance of accounting for the bend- 
ing moment in the joint due not only to the eccentricity of the loading but also 
to the deformation of the adherends as the load is applied. This interdepen- 
dence of the bending moment, load, ard deformation makes the problon geometri- 
cally nonlinear. Since the SPAR program is linear, an approximation to the 
nonlinear behavior is computed using SPAR by first obtaining a linear solution 
with the bending moment based on the undeformed shape. The initial stress- 
stiffness matrix associated with the applied load and the resulting deformations 
is calculated and then added to the initial stiffness matrix. Next, the linear 
SPAR analysis is rerun with the modified stiffness matrix to give a one-step 
approximation to the nonlinear problem. This process is repeated until conver- 
gence of the solution is achieved. It was found that only one recomputation was 
required to obtain convergence for the lap-joint analysis under consideration. 


Nonlinear Finite-Difference Analysis 

In order to verify the accuracy of the nonlinear finite-element approach, 
the single lap joint was analyzed by using the nonlinear finite-difference code 
STAGS described in reference 27. STA®I was developed primarily for the non- 
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linear analysis of shell- or plate structures and includes geometric as well 
as material nonlinearity. 

The finite-difference mesh used to analyze the single lap joint is shown 
in figure 9 (b) . Analyses using different mesh sizes indicated that a total of 
1019 finite-difference stations were adequate to give converged results. 
Finite-difference nodes were concentrated in the critical areas where large 
stress gradients occur, such as at the end of the adhesive and in the surrovmd- 
ing adherends. (See enlarged detail in fig. 9(b).) 


Comparison of Numerical Results 

The results of the finite- element and finite-difference analyses for non- 
dimensional shear and peel stress distributions along the center line of the 
adhesive are shown in figures 10 and 11. Both linear and nonlinear results 
are presented. The solid and dashed curves were obtained with the finite- 
element and finite-difference analyses# respectively. The two analyses agree 
well, especially in the critical region near the free edge of the adhesive 
shown by the expanded scale. 

A comparison between the linear and nonlinear results shows that for a 
single lap joint the change in bending moment due to the lateral deflection of 
adherends should be taken into account. The maximum shear stresses along the 
center of the bond line in figure 10 differed by 28 percent, and the maximum 
peel stresses in figure 11 differed by 35 percent. 


COMPARISON OP GQLAND-ifflISSNER AND FINITE- ELEMENT RESULTS 

In figures 12 and 13, the shear and peel stresses for 175 kN/m 
(1000 Ibf/in.) applied load obtained with the Goland-Reissner analysis taken 
from figures 4 and 5 are compared with the stress distributions along the mid- 
plane of the adhesive computed using both the linear and nonlinear finite- 
element procedures. Because of the high stiffness of the adhesive used in the 
current study, the lap-joint configuration under consideration is outside the 
suggested range of applicablity of the GolaM-Reissner analysis as stated in 
reference 4 and slightly outside the relaxed range of applicability suggested in 
reference 28. However, with the exception of a small region near the overlap 
edge, agreement is excellent. The Goland-Reissner analysis is quite accurate in 
predicting the shear stress distribution at the center line of the adhesive over 
at least 95 percent of the bond line. 

As seen in figure 12, the maximum shear stress concentrations computed by 
the nonlinear analyses differ by 12 percent, with the Goland-Reissner analysis 
giving the higher value. Both of these maximum values are considerably below 
the maximum shear stresses predicted by the corresponding linear analyses. The 
free-edge effect is shown in the figure and is restricted to a region within 
2 percent of the overlap length. 

As seen in figure 13, the maximum peel stress concentrations predicted by 
the Goland-Reissner and finite- element nonlinear analyses differ by approxi- 
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mately 25 percent/ the Goland-Reissner analysis giving the higher value. The 
Goland-Reissner analysis is less accurate throughout a lar^r part of the over- 
lap region in the prediction of the peel stress distribution than for the pre- 
diction of the shear stress distribution, but still gives good results over 
approximately 70 percent of the bond line. 

As seen from these comparisons, the lack of satisfaction of the adhesive 
free-edge boundary condition in the Goland-Reissner analysis has a minor effect 
in the prediction of the adhesive stress distributions along the center line 
of the adhesive except near the overlap edge. The form of the stress distribu- 
tions and the prediction of the stress level are sufficiently accurate for 
determination of the qualitative behavior of single lap joints. Thus, the 
stress equations (A8) and (AW) provide an efficient procedure for a qualitative 
evaluation of the influence of different geometries and material properties on 
the center-line adhesive stresses. 

The Goland-Reissner analysis assumes that the adhesive shear and peel 
stresses are constant through the thickness and neglects the in-plane direct 
stress. The shear, peel, and in-plane direct stress distributions at three 
locations through the adhesive thickness are plotted along the bond line in 
figures 14 to 16 and are compared with the results from the Goland-Reissner 
analysis. Over most of the bond line the shear and peel stresses show very 
little variation through the thickness. However, there is considerable varia- 
tion through the thickness near the free edge of the adhesive due to the 
influence of the free-edge boundary condition and proximity to the reentrant 
corner. The direct stress in the adhesive is of the order of 10 percent of 
the applied adherend normal stress along most of the overlap but increases 
considerably in the vicinity of the free edge. 


LAP-JOINT DETAILED STRESS DISTRIBUTIONS 

Detailed stress distributions presented in the following sections can be 
used for comparison with other solution techniques and analyses. Stress dis- 
tributions are presented through the thickness and along the length of both the 
adhesive and the adherend in figures 14 to 23. Since the finite-element analy- 
sis is an elastic analysis, a stress singularity exists at the reentrant corner 
where the free edge of the adhesive joins the adherend. The localized nature 
of the singular point is most evident in the stress distributions of the adhe- 
sive in figures 15 and 18 and in the adherend in figures 22 and 23. A further 
refinement of the finite-element modeling of the adherend and adhesive in the 
vicinity of the reentrant corner was made; and, althoigh the maximum computed 
stress values increased rapidly ais the singular point was approached, the vedues 
presented in the figures are at a sufficient distance from the singularity to 
remain unchanged. Although the stress data are presented in nondimensional form, 
it must be realized that the stresses vary nonlinearly with applied load ard were 
calculated using a tensile load T of 175 kN/m (1000 Ibf/in.) applied to a lap 
joint with dimensions given in figure 1. Use of nondimensional stress values 
for other load levels or other specimens would not be correct, as can be seen 
from figures 4, 5, 7, and 8. 
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Nondimens ional stress distributions in the adhesive are presented in fig- 
ures 14 to 19. Shear, peel, and in- plane stress distributions along the length 
of the adhesive at different positions through the thickness are shown in fig- 
ures 14, 15, and 16, respectively. Shear, peel, and inrplane stress distribu- 
tions across the adhesive thickness at several different positions along the 
joint length are shovm in figures 17, 18, and 19, respectively. 

Nondimensional transverse and axial stress distributions in the adherend 
are shown in figures 20 to 23. Figures 20 and 21 show the transverse and axial 
stress distributions through the adherend thickness at various positions along 
the adherend, and figures 22 and 23 show the stress distributions along the 
adherend at various positions through the thickness. The transverse normal 
stresses (figs. 20 and 22} in the adherend are near zero except near the free 
edge of the adhesive where large values and variations occur both along and 
across the adherend. The large transverse tensile stresses and gradients prove 
especially detrimental for filamentary composite single lap joints due to the 
low interlaminar strength of most ccmposlte laminates. 

The axial stresses vary linearly through the adherend thickness (fig. 21 ) 
except at the edge of the overlap. Near the center of the joint (x/il = 0, 
fig. 23) , the axial stresses are approximately constant through the adherend 
thickness and are approximately 50 percent of the average applied stresses. 

Near the applied load, the axial-stress variation through the adherend thick- 
ness is small, and the stress values are approaching the average applied 
stresses. Along the adherend, the variations in the axial stress through the 
adherend thickness become larger near the overlap edge, thus reflecting the 
distribution in bending-moment trend shown in figure 3. 


CRITIQUE OF ASTM D 1002-72 ADHESIVE SHEAR-STRENGrTH TEST 

since the shear strength of an adhesive must be known before a joint design 
can be made, a practical inexpensive test procedure to evaluate the shear 
strength is desirable. The analyses performed in this paper are concerned with 
a lap-joint configuration similar to that suggested by ASTM D 1002-72, which is 
the most popular procedure for the evaluation of the shear strength of adhesive 
systems. The results suggest, in agreement with findings of others (refs. 19 
and 21), that this test is inadeguate even in a qualitative sense. A test 
specimen using an adhesive weak in peel strength but with high shear strength 
could fail at a load below that required to fail a test specimen bonded with an 
adhesive strong in peel strength but weak in shear strength. Thus, the evalua- 
tor could be led to an erroneous conclusion as to the relative shear strengths 
of the two adhesives. It appears that other established test techniques 
(ASTM E 229-70 or the thick- adherend test given in ref. 29), although lacking 
the simplicity and economy of the ASTM D 1002-72 test, should be used for evalu- 
ating adhesives or consideration should be given to the development of a simple 
economical test specimen which can replace the ASTM D 1002-72 specimen. 
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CONCLUDING REMARKS 


A linear elastic finite- element analysis of a bonded single lap joint is 
modified to account for the nonlinear behavior of the joint due to the eccen- 
tricity of the load path. The influence of the nonlinear behavior on the 
stress distribution in the adhesive and adherend is sufficiently large so that 
a sizable error in the prediction of the adhesive shear and peel stress dis- 
tributions can occur if the nonlinear effects are not included. This nonlinear 
behavior was shown experimentally in a photoelastic study of the single-lap- 
joint adherend. 

Results from the classic closed-form approximate analysis of Goland and 
Reissner for an elastic single lap joint/ which does take into account the 
flexibility of the adherends and joint rotation, are compared with the finite- 
element solution for a lap-joint configuration similar to the ASTM D 1002-72 
standard test specimen for evaluation of the shear strength of adhesives. 
Results from this comparison indicate that the Goland-Reissner solution is suf- 
ficiently accurate to be used in the qualitative prediction of the influence 
of parametric variations on the shear and peel stresses along the center of 
the adhesive bond line. The closed-form solution does not satisfy the adhesive 
free-edge shear boundary condition, but the influence of the free edge is con- 
fined to a small zone and has little influence on the computed maximum value 
of the shear or peel stress occurring along the center of the adhesive bond 
line. Inspection of the adhesive stress distribution obtained with a finite- 
element analysis shows that the assumption of constant shear ard peel stress 
through the adhesive thickness used in the Goland-Reissner solution is valid 
over the interior of the lap region but deteriorates in the vicinity of the 
free-edge reentrant corner where a two-dimensional elastic analysis predicts 
singular stresses. 

The method of evaluating the interned rncment in the adherend used in the 
Goland-Reissner approach, which assumes that the joint deforms as a plate in 
cylindrical bending, is shown by both the finite-element results and a photo- 
elastic study to be accurate even near the region of the overlap for an adhe- 
sive with a modulus much less than the modulus of the adherend. 

Detailed stress distributions at several locations in both the adherend 
and adhesive are presented to provide a basis for comparison with other solu- 
tion techniques. Inspection of the stress distribution in the ASTM D 1002-72 
lap-joint configuration shows that this test specimen is inadequate for deter- 
mining the shear strength of adhesives even in a qualitative sense. 


Langley Research Center 

National Aeronautics and Space Adainistration 
Hampton, VA 23665 
July 27, 1979 



APPiMDIX A 


DERIVATICai OF GOLAND-KEISSNER EQCJATIONS FOR ADHESIVE SHEAR 

AND PEEL STRESS 

The development of the shear and peel stress equations is similar to that 
of reference 4 , but with plate behavior being assumed consistent throughout the 
development and, in addition, small terms formerly Ignored being retained for 
ccanpleteness. Wherever terms differ from corresponding terms in reference 4, 
these terms are underlined with a dashed line. 

The sign convention used throughout this analysis is shewn in figure 24. 

A differential element of a slice through the bonded region is shown in fig- 
ure 25 where the adhesive is assumed to act as a spring with both shear and 
transverse extensional stiffness so that stress is assumed constant through the 
thickness. The integrated internal stresses (stress resultants) are applied 
on the vertical cut, thiis reducing the analysis to a one-dimensional plane- 
strain analysis. 

From equilibrium considerations. 



(Ala) 


(Alb) 


(Ale) 


(Aid) 


(Ale) 


(Alf) 
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APPENDIX A 


where the subscripts D and L refer to quantities in the upper and lower 
adherends of figure 25. 

The assumption of plane strain leads to the load- displacement relation- 
ships of cylindrically bent plate theory 


d^Wy 

dx2 


Mu 

D 


d2wL Ml 

dx2 ^ 


where the flexural stiffness of the adherend is given by 


Et^ 

D = 

12(1 - v2) 


(A2a) 


(A2b) 


The reduction of the flexural stiffness due to the existence of the adhe- 
sive bond line is ignored, and shear deformations in the adherend are neglected 
as being small compared with flexural deformations. The longitudinal strain 
at the adhesive-adherend interface due to flexure is given by 



t % 
2 D 



so that the total strain in terms of forces and moments becomes 

duu Tu t Mu 

dx A 2D 

duL Tl t Ml 

dx A 2D 


(A3a) 


(A3b) 


13 



APPENDIX A 


where 



Since the adhesive shear and transverse normal stresses are assumed con- 
stant through the adhesive thickness, then 



(A4a) 



(A4b) 


If only the overlap region is considered, the stress resultant continuity condi- 
tions can be given as 


Vl = Ml = Tl = 0 

Vu = Vo; Mu = Mq; Tu = To 


Vu = Mu = Tu = 0 

Vl = Vo; Ml = -Mo; Tl = To 


(X = -c) 


(X = c) 


(A5a) 


(A5b) 


where quantities with subscript o are the internal stress resultants applied 
by the adjoining adhe rends. 

Differentiating equation (A4a) twice and equations (A3) once, combining 
the resultant equations to remove Ul and Uu, and substituting for the deriva- 
tives of M and T by using equations (Ala) to (Aid) result in the following 
equation; 


d^T 

dx^ 
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APPENDIX A 


On differentiating this equation once and using equations (Ale) and (Alf) to 
remove V, the ecpiation governing the distribution of shear stress in the adhe- 
sive is found to be 


dT 

<3x3 Vc/ dx 


(A6) 


where 


62 = 


8Ga(l - v2) 


HEt 



and 




n 

t 


The equation is subject to the following boundary conditions found from equa- 
tions (A3) to (A5) and equilibrium of an adherend in the lap region: 


dT 

dx 


(-C) 


-Ga(l - / 

— — —To 

TlEt \ 



(A7a) 



Gad - v2) 
riEt 



+ 



(A7b) 


Additionally from equilibrium of external forces, 

+c 

T dx = To {A7c) 

c 
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APPENDIX A 


The solution of the system of equations (A6) and (A7) is 



4 + 3C 


6 


(1 + 3 k) 


sinh 5 


cosh(~j + 3(1 - k +_?) 


(A8) 


where 


■^av 


- 

2c 


and 


2Mq 
~ ^ 


which reduces to the solution of reference 4 if ? = 0 and the term 1 - v2 is 
removed from the equation for 5 2 , 

Differentiating equation (A4b) three timeS/ combining with equations (A2) 
differentiated once to remove w, and combining with equations (Ale) and (Aid) 
result in the following: 


d^a 

dx3 


Ea 

= — (Vu - Vl) 
Dh 


On combining this equation with equations (Ale) and (A1f), the eqpiation govern- 
ing the peel stress in the adhesive is found to be 


d^a 2Ea 

— + — a a 0 (A9) 

dx4 Dn^ 


Equation (A9) is subject to the following boundary conditions found from equa- 
tions (A2) , (A4) , Snd (A7)-: 


16 



APPENDIX A 


d^a 

dx2 


-c 


d2a 

dx2 


E, 


no 




d^a 

dx2 


-c 


d3a 

dx2 


Ea 

= — 


Ic 


r\D 


The solution of the system of equations (A9) and (AT 0) is 


a 


't\^ 1 


^app \cy ^ 


\ /Xx\ /Xx\ 
+ k'A sinh A sin A j sinhl — 1 sinj— J 


/r2^^k \ /Xx\ /Xx\ 

+ I — ^ — + k ' X cosh X cos X I cosh( — | cosj - 


\c 


\c 


where 


app 


X = c4, 


2r)D 


A = -(sin 2X + sinh 2X) 
2 


r-| ® cosh X sin X + sinh X cos X 


r 2 = sinh X cos X - cosh X sin X 


CVo 


k' = 


tTr 


(AT Oa) 


(A1 Ob) 


(All) 
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The values of Tq, Mq, and Vq must be evaluated to complete the solu- 
tion. To obtain approximations of these internal stress resultants, the lap 
joint is assumed to behave as a wide beam loaded in in-plane tension with a 
discontinuous neutral surface, as shown in figure 26. This assumption is of 
questionable accuracy near the lap region where its accuracy depends to a large 
degree on the relative flexibility of the adhesive. To examine this assumption, 
two photoelastic lap-joint models with stiff and flexible adhesives were ana- 
lyzed. The results of this study are given in appendix B. 

The cross-sectional moments in regions 1 and 2 can be given by 


M-| = T cos <t> (x-| tan - w-| ) 


M 2 = T cos <j> [j[X 2 - c) tan <j> - W 2 ^ 

so that the beam equations governing the total lap-joint deformation can be 
given as 


<3^1 T cos (|> 


dx^ 


(Xi tan <1) - W] ) 


(0 < xi < h (A12a) 


d^w 

dx2 


2 T cos 4> 
8D 


[<X 2 - c) tan <l> - W 2 ^ (0 < X 2 c) (AT 2b) 


subject to the boundary and continuity conditions 


"i l o “ ° 

«2lc “ 0 

U = ''2lo 


dw-] 

dx 


dw2 

dx 


(A1 3a) 
(AT 3b) 
(A13c) 

(A1 3d) 


The presence of terms containing the product of T and w make equations (A12) 
nonlinear with respect to the applied load. In addition, from geometric consid- 
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erations and simple plate theory, the stress resultants at the interface 
the lap region and the single adherend can be given by; 


Tq = T cos (j) 


1 d^w-j 


Mo = - - 


^ dx2 


Vo = 


1 d^wi 

D dx3 


Solution of the system of equations (AT 2) to (A14) yields 


ft + n\ cosh(U2C) sinh<uiS-) 

Mq = T cos <1) ' ' • ■ . - - — 

\ 2 / U2 sinh(uiX) cosh(U2C) + U] cosh{ui5«) sinh(U2c) 


Vo = T| 


rt + h) 


cos <t> 


uiU 2 cosh(U 2 C) cosh (ui S') 


U 2 sinh(uiA) cosh(U 2 C) + u-] cosh(u-|ii) Sinh(U 2 C) 


where 


cos <|> = 


a + c 


(Z + e) 2 + 


ft + 


ui = 


|T cos 


U2 = 


T cos (j) 
8D 


between 


(A14a) 
(AT 4b) 


(AT 4c) 


(AT 5a) 


(AT 5b) 
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The solution of the stress resultants reduces to the values of Tq, Mq, and 
Vq of reference 4 for cos <j) = 1 . 

If the adherends are assumed rigid so that D approaches infinity, the 
applied moment and shear distributions in the adherends reduce to 


Ml 



(A1 6a) 


Vi 



cos <j> 


(A1 6b) 


and the adhesive shear and peel stresses given by equations (a 8) and (All) are 
linearly related to the applied tensile-force resultant. 


/ 
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PHOTOELASTIC EXAMINATI(»I OF LAP-JOINT ADHEBENDS 

Three photoelastic models shown in figure 27 were examined to evaluate 
the moment distribution in the adherend in the region outside the overlap for 
single lap joints with either stiff or flexible adhesives. The adherends of 
the models were fabricated from an epoxy material (PSM-5) . The models were 
sized to approximate the geanetric ratios of the ASTM D 1 002-72 test specimen. 
Several photoelastic studies have been made previously of lap joints (refs. 20 
and 30 to 32 ) , but in each case the models considered had impractical geometri- 
cal ratios. The first model, shown in figure 28(a), is an integral model repre- 
sentative of a lap joint with an adhesive of the same material as the adherend. 
The second lap-joint model (fig. 28(b)) is an actual joint bonded with PMG-1 
adhesive, which has properties similar to the adherend and cures at room tem- 
perature. This model has a stress-distribution pattern \^ry similar to the 
integral model. The third model (fig. 28(c)) is a lap joint with a soft epoxy 
material (PS-3) bonded between the adherends to represent an adhesive with shear 
and extensional moduli well below the adherend moduli. The material properties 
and required stress optic constant for the PSM-5 photoelastic material are given 
in table II. All models were loaded to 133 N (30 Ibf). 

The experimental mcments are determined by assuming that the adherend 
deforms as a simple beam subject to in-line tension so that the principal 
stresses determined frcm the isochromatic fringes are assumed directed along 
the axis of the adherend. Inspection of the isoclinics showed this assumption 
to be nearly correct everywhere in the adherend outside the overlap region for 
the model with the flexible adhesive and everywhere except close to the overlap 
in the integral and stiff adhesive models. The axial stresses are determined 
at several points across the adherend model width for each axial location. A 
best-fit straight line is computed using a least-squares procedure and the stress 
determined at one boundary by extrapolation. The results for three locations 
along the adherend in the integral model are shown in figure 29. The model with 
the stiff adhesive exhibited the same behavior as the integral model, as seen 
by comparison of the isochromatics shown in figures 28(a) and 28(b). The stress 
distributions are nearly linear across the width of the model at each location 
along the adherend except near the overlap region. For the model with the soft 
adhesive, the stress distributions are nearly linear at all locations outside 
the overlap region. The moment is computed from the approximate beam-rod 
equation 



where I is the moment of inertia of the cross section about the cross- 
sectional centroid, a is the area of the adherend cross section, Ojj is the 
extrapolated stress at the boundary, t is the width of the adherend (fig. 29) , 
and P is the applied load. The evcduated mcanent distribution in the adherend 
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of each model near the overlap region is shown in figure 30 and is compared with 
the moment distribution calculated frcm the Goland-Reissner analysis for the 
stiff emd flexible adhesive models. 

The Goland-Reissner analysis neglects the adhesive in the ccmputation of 
the monent and shear-force resultants in the adherend. As shown in figure 30, 
the agreement between the canputed and experimentally determined moments for 
the flexible adhesive model is good everywhere in the adherend (within 9 per- 
cent) . Thus, the error due to neglect erf the bending stiffness of the adhesive 
in the computation of the mcanent is slight. The agreement is also good for 
the stiff adhesive models, except close to the lap edge where large transverse 
normal stresses and shear stresses begin to dominate - resulting in 15 percent 
error. Neglect of these stresses in the computation of the moment in the vicin- 
ity of the lap can lead to same error. Further, and perhaps more significant, 
a comparison of the experimental results with the moment calculated assuming 
a rigid adherend, which is the assumption tacitly made using a linear analysis, 
shows that the difference in moment due to the flexibility of the adherend is 
considerable and predicted well by the Goland-Reissner analysis. This compari- 
son demonstrates experimentally the geometric nonlinear effect in the calcula- 
tion of the adherend moment at the edge of the lap Mq. 
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TABLE I.- LAP-JOINT MATERIAL PROPERTIES 


Adherend 

Young' s modulus, GPa (psi) . 110 (16 x 10®) 

Shear modulus, GPa (psi) . . . . .... . 43. (6.2 x 10®) 
Poisson's ratio ............. 0.3 


Adhesive 
17.9 (2.6 X 10®) 
3.2 (0.46 X 10®) 
0.19 


TABLE II.- PHOTOELASTie MATERIAL PROPERTIES 



PSM-5 

PMC-1 

PS-3 


Young' s modulus. 

3.1 (450 X 103) 

2.9 (425 X 103) 


103) 

GPa (psi) 

0.21 (30 X 

Poisson's ratio . . . . . 

0.36 

0. 36 


0.3 

Fringe value, kPa-m 





per fringe (psi-in. 
per fringe) ...... 

10 (57.4) 
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13.97 ( 5 . 50 ) 



Pigtire l.~ Single-lap-joint gecanetry. Dimensions are given in 

centimeters (inches). 




175kN/m 



Figure 3.- Moment distribution in adherend. 




Figure 4.- Adhesive shear stress distribution at several load levels. Goland-Reissner 

analysis; c = 0.64 cm (0.25 in.); = T/2c. 















(a) SPAR (ref. 25) 


finite- element mesh. 



(b) STAGS (ref. 27) finite-difference mesh. 

Figure 9.- Pinite-elonent and finite-difference model of single lap joint. 
Dimensions are given in centimeters (inches). 






Figure 11.- Conparison of peel stress distributions along center line 
of adhesive given by finite-difference and finite-elenent analyses 


T = 175 kN/m (1000 Ibf/in.); a^pp = T/t. 



0 2 
^ See expanded view * 


Figure 12.- Comparison of shear stress distributions along center line of 
adhesive given by Goland-Reissner analysis and finite-element analysis 
T = 175 kN/m (1000 Ibf/in.); = T/2c. 
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Figure 15.- Peel stress distributions at various positions through adhesive. 
T = 175 kN/m (1000 Ibf/in.); Ogpp = T/t. 
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Figure 16.- In-plane stress distributions at various positions through 
adhesive. T = 175 kN/m (1000 Ibf/in. ); <^app “ T/t. 
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Figure 25.- Forces on differential element 
shovm in positive direction. 
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Figure 28 
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,“ Concluded, 
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Transverse distance through adherend, y/t 

Figure 29.- Photoelastic measurement of adherend stress distribution near 
overlap. Stiff adhesive model; CTapp “ P ^ 33 N (30 Ibf ) ; 

a = 0.32 cm^ (0.05 in^) ; H = 22.0 cm (8.65 in.); c = 4.32 cm 
(1.70 in.); t = 0.86 cm (0.34 in.). 
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Figure 30.- Mcanent distribution along adherend near overlap 

T = 1 33 N (30 Ibf ) . 



